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' Abstract. We review the major progress of the past decade concerning our 

understanding of the nucleon-nucleon interaction. The focus is on the low-energy 
' region (below pion production threshold), but a brief outlook towards higher 

energies is also given. The items discussed include charge-dependence, the precise 
value of the nNN coupling constant, phase shift analysis and high-precision NN 
1^"^ ' data and potentials. We also address the issue of a proper theory of nuclear forces. 

' Finally, we summarize the essential open questions that future research should be 
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1. Introduction 



The nuclear force has been at the heart of nuclear physics ever since the field was born 
in 1932 with the discovery of the neutron by Chadwick [1]. In fact, during the first few 
decades of nuclear physics, the term 'nuclear forces' was often used as synonymous 
for nuclear physics as a whole [2] . There are good reasons why the nuclear force plays 
^ I such an outstanding role. 

The interaction between two nucleons is basic for all of nuclear physics. The 
traditional goal of nuclear physics is to understand the properties of atomic nuclei in 
' terms of the 'bare' interaction between pairs of nucleons. With the onset of quantum- 

H I chromodynamics (QCD), it became clear that the nucleon-nucleon (NN) interaction 

is not fundamental. Nevertheless, even today, in any first approach towards a nuclear 
structure problem, one assumes the nucleons to be elementary particles. The failure or 
success of this approach may then teach us something about the relevance of subnuclear 
degrees of freedom. 

The NN interaction has been investigated by a large number of physicists all 
over the world for the past 70 years. It is the empirically best known piece of strong 
interactions; in fact, for no other sample of the strong force a comparable amount of 
experimental data has been accumulated. 

The oldest attempt to explain the nature of the nuclear force is due to Yukawa [3] . 
According to his theory, massive bosons (mesons) mediate the interaction between two 
nucleons. This idea spawned the sister discipline of particle physics. Although, in the 
light of QCD, meson theory is not perceived as fundamental anymore, the meson 
exchange concept continues to represent the best working model for a quantitative 
nucleon-nucleon potential. 
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Historically, it turned out to be a formidable task to describe the nuclear force 
just phenomenologically, and it took a quarter century to come up with the first semi- 
quantitative model [4] — in 1957. Ever since, there has been substantial progress in 
experiment and theory of the nuclear force. Most basic questions were settled in the 
1960's and 70's such that in recent years we could concentrate on the subtleties of this 
peculiar force. 

In this topical review, we will report the chief progress of the past decade. The 
focus will be on the low-energy region (below pion production threshold) . Summaries 
of earlier periods and a pedagogical introduction into the field can be found in 
references [5, 6]. In the 1990's, major issues concerning the NN interaction have 
been: 

• charge-dependence, 

• the precise value of the ttNN coupling constant, 

• improved phase shift analysis, 

• high-precision NN data, 

• high-precision NN potentials, 

• QCD and the nuclear force, 

• NN scattering at intermediate and high energies. 
We will now review these topics one by one. 

2. Charge dependence 

By definition, charge independence is invariance under any rotation in isospin space. A 
violation of this symmetry is referred to as charge dependence or charge independence 
breaking (CIB). Charge symmetry is invariance under a rotation by 180'^ about the 
y-axis in isospin space if the positive z-direction is associated with the positive 
charge. The violation of this symmetry is known as charge symmetry breaking (CSB). 
Obviously, CSB is a special case of charge dependence. 

CIB of the strong NN interaction means that, in the isospin T = 1 state, the 
proton-proton (T^ = +1), neutron-proton (T^ = 0), or neutron- neutron {T^ = — 1) 
interactions arc (slightly) different, after electromagnetic effects have been removed. 
CSB of the NN interaction refers to a difference between proton-proton (pp) and 
neutron- neutron (nn) interactions, only. The charge dependence of the NN interaction 
is subtle, but in the ^Sq state it is well established. The observation of small charge- 
dependent effects in this state is possible because the scattering length of an almost 
bound state acts like a powerful magnifying glass on the interaction. 

The current understanding is that — on a fundamental level — the charge 
dependence of nuclear forces is due to a difference between the up and down quark 
masses and electromagnetic interactions among the quarks. A consequence of this are 
mass differences between hadrons of the same isospin multiplet and meson mixing. 
Therefore, if CIB is calculated based upon hadronic models, the mass differences 
between hadrons of the same isospin multiplet, meson mixing, and irreducible meson- 
photon exchanges are considered as major causes. For reviews on charge dependence, 
see references [7, 8, 9]. We will now summarize recent developments (that are not 
contained in any of these reviews). 
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2.1. Charge symmetry breaking 

2.1.1. Experiment. As discussed, the scattering Icngtlis in the state for pp, np, 
and nn scattering (denoted by app, a„p, and a„„, respectively) are the best evidence 
for the charge-dependence of nuclear forces. While we have well-established values for 
app and Unp since many decades, the neutron-neutron scattering length continues to 
be a tough problem. The basic reason for this is that, so far, we have not been able 
to conduct any direct measurements of a„„ using free neutron- neutron collisions [10] . 
All current values are extracted from multi-particle reactions the analyses of which 
are beset with large theoretical uncertainties. The processes that are believed to have 



the smallest uncertainties are 

li~+d^Vfi + n + n, (1) 

n'+d^j + n + n, (2) 

n + d ^ p + n + n. (3) 



While, there are no data on the first reaction, the other two processes have been 
studied repeatedly. In 1998, a very carefuU study of the 7r~ induced reaction was 
published [11] and, in 1999, a renewed thorough investigation of the neutron induced 
process was accomplished [12], yielding results that are in perfect agreement, namely. 



D{n-,n'y)n [11] : a„„ = -18.50 ± 0.53 fm, (4) 
D{n, nnp) [12] : a„„ = - 18.7 ± 0.6 fm, (5) 

which can be summarised by 

a„„ = -18.6±0.4fm. (6) 
Correcting for the neutron-neutron magnetic interaction, the pure nuclear value is: 

a^„ -18.9±0.4fm. (7) 



This summarizes the status by the end of 1999. Unfortunately, this is not the happy 
end of the story that everybody had hoped for. To properly discuss the new (and old) 
problems, we will first provide more details concerning the two types of reactions for 
which experiments have been conducted. 

Over the past 20 years, there have been three independent studies of the reaction 
TT^ + 7 + n + n. In one case [13], only the 7 spectrum was measured, while 

in the other two cases [14, 11], kinematically complete experiments were performed 
measuring the 7 and a neutron in the final state. The results are: 

a„„ = - 18.60 ± 0.34 (stat.) ± 0.26 (syst.) ± 0.30 (theor.) fm 

= - 18.60 ±0.52 fm [13], (8) 

a„„ = - 18.70 ± 0.42 (stat.) ± 0.39 (syst.) ± 0.30 (theor.) fm 

= - 18.70 ±0.65 fm [14], (9) 

a„„ = - 18.50 ± 0.05 (stat.) ± 0.44 (syst.) ± 0.30 (theor.) fm 

= - 18.50 ±0.53 fm [11]. (10) 

Owing to the high spatial resolution of the gamma ray detector in reference [11], it 
was possible to assess the systematic errors due to uncertainties in the modelling of 
the stopped pion distribution in the target and in target vertex reconstruction in the 
Monte Carlo simulation. Therefore, the systematic uncertainties of the kinematically 
complete studies are now much better understood, and a very high statistical accuracy 



Nucleon-nucleon interaction 



4 



in reference [11] makes the experimental uncertainty comparable to the theoretical one 

in the extraction of a„„ from the reaction D{Tr~ ,jn)n. The combined result from all 
three studies gives the new world average for the D(7r^,7)nn reaction 

a„„ = - 18.59 ± 0.27 (exper.) ± 0.30 (theor.) fm 

= - 18.59 ±0.40 fm. (11) 

In summary, the reaction tt" + d^j + n + n appears to be in good shape. 

Unfortunately, we cannot say the same about the neutron-induced deuteron 
breakup process. Until the recent investigation by the TUNL group, Gonzalez-Trotter 
et al. [12], all studies of the reaction n + d — > p + n + n gave for a„„ values that 
differed from that obtained from the Din^ ,jn)n process. Theoretical uncertainties 
in extracting a„„ from the neutron- induced deuteron breakup are much larger, as we 
will explain now. 

First, in reactions with more than two nuclcons in the final state, three miclcon 
forces (3NF) modify the cross section. It was suggested [15] that the 3NF is the 
reason why a„„ extracted from the D(n, nnp) process differs from that obtained from 
the D(7r^,7n)n reaction. The 3NF arc a natural consequence of strong interactions. 
Therefore, 3NF do exist, but the question is how significant they are and, in particular, 
do they affect a specific configuration of the nd breakup wherefrom one extracts Unn ■ 
There are several indications for possible 3NF effects in nuclear physics: •'ll binding 
energy, nuclear matter binding energy, ^He binding energy and negative parity excited 
states, ^He and ^He one-body density distributions, binding energies and radii of some 
nuclei, ^''O magnetic moment form factor, nd capture, Ay in elastic nd scattering, and 
space star, final state interaction (FSI) and quasifree scattering (QFS) configurations 
in the nd breakup [7, 16, 17]; but none of them provided conclusive information on 
3NF. It was possible to reconcile all values of a„„ extracted from the studies of the 
reaction D{n,nnp) in the energy domain of 10 to 50 MeV with those obtained from 
the D(7r~,7n)n process using the Fujita-Miyazawa 3NF [7]. However, the reanalyses 
of these nd breakup processes [18] gave values that differed considerably from those 
quoted by the original authors. Though any reanalysis is clouded by the lack of all 
relevant information, the main reason is the fact that original analyses used simple S- 
wave separable NN potentials, while the re- analyses were done using the rigorous three 
body theory of Glockle et al [16]. Obviously, the claimed theoretical uncertainties in 
the original papers were underestimated. 

Second, the magnetic interaction modifies the value of the ^Sq scattering length 
extracted from the neutron induced deuteron breakup. It was shown [19] that for 
the neutron-pickup configuration in the neutron-induced deuteron breakup leading to 
the nn FSI there is a magnetic interaction in the ^ 5*0 state which is repulsive thereby 
decreasing the absolute value of a^n- Depending on the NN potential (hard core or soft 
core), impulse approximation estimates of the effect of the magnetic interaction in the 
pickup configuration changes a„n from —18.5 fm to —17.2 or even —16.4 fm [19]. The 
correction for the knockout configuration has the opposite sign since it is dominated 
by the magnetic interaction between a neutron and a proton in the ^Sq state. The 
situation is more complex for the neutron-proton FSI, since the np FSI occurs in the 

and "^^i states. 

The determination of a„„ by Gonzales-Trotter et al [12] has two characteristic 
features: first, it uses the rigorous theory [16] including, in addition to several realistic 
NN potentials, also the Tucson-Melbourne 3NF, and second, it performs a high- 
accuracy comparison of neutron- neutron and neutron-proton FSI in the ^Sq state by 
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measuring cross sections of the reaction D{n,nnp) for identical kinematic conditions 
(the angle of the two emitted nucleons interacting in the ^ 5*0 final state is 28 to 43 deg) 
at the incident neutron energy of 13 MeV. Therefore, the neutron-proton scattering 
length, a„p, becomes the standard for determining o„„. By comparing the extracted 
value for a„p and its uncertainty, it was possible to set an upper limit of 0.2 ± 0.6 fm 
on any possible effects due to 3NF influencing the extracted value of a„„. Of course, it 
is possible that the effect of the magnetic interaction discussed by Slobodrian et al [19] 
arc negligible in the energy/angular region studied by Gonzalcz-Trottcr et al [12]. On 
the other hand, it should be stressed that the rigorous calculations by Glockle et 
al [16] do not include the electromagnetic interaction, and that there are now many 
indications of the shortcomings of the Tucson-Melbourne 3NF. 

In the year of 2000, a new study of the neutron-neutron FSI in the D{n,nnp) 
reaction at the incident neutron energy of 25.3 MeV was published by the neutron 
group at Bonn [20]. The data were analyzed using the rigorous theory [16]. The 
extracted value is. 



which is in drastic disagreement with the result of the TUNL group [12] published in 
1999 and also with those obtained from the D{-!t~ ,^n)n process. 

While most of the previous kinematically complete studies of the reaction 
D(n, nnp) employ a thick, active deuterated target measuring the energy of the proton, 
and detecting two neutrons at nearly the same angle on the same side of the incident 
neutron, this recent measurement [20] uses a thin deuterium target and detects a 
neutron at 9„ = —55.5 deg and a proton at Qp = 41.15 deg. The advantage of this 
geometry is the reduction of the strong cross talk between neutron detectors and the 
reduction in losses from neutron multiple scattering. This geometrical configuration 
has the added advantage that the locus contains np QFS besides nn FSI and, therefore, 
provides a built-in normalization. Indeed, normalizing the data to np QFS yields a 
very similar value: 



Neither the use of different NN potentials nor the inclusion of the Tucson-Melbourne 
3NF in the rigorous calculation produces noticeably different results for a„„. This 
geometry at this incident energy is the region where the 3NF effect of the Tucson- 
Melbourne potential is very small. The preliminary result by the Bonn group — using 
the same incident energy of 25.3 MeV — gave a good fit to the np FSI spectrum using 
dnp = — 24 fm. 

The disagreement between the two most recent studies, Gonzalez- Trotter et al 
(TUNL) [12] and Huhn et al (Bonn) [20], opens again the problem of how completely 
do we understand the interactions involved in the three nucleon problem, specifically 
the 3NF. It also suggests that additional experimental studies at different incident 
energies and at different angles might be useful in resolving the problem. 

When we use for a,,,,, the value obtained from the D{'K^,"f)nn studies 
[equation (11)], correct it for the magnetic moment interaction [equation (7)], and 
compare it to the corresponding pp value [8]: 



a„„ = -16. 27 ±0.4 fm. 



(12) 



a„„ = -16.06 ±0.35 fm. 



(13) 



a^p = -17.3 ±0.4 fm, 
then charge-symmetry is broken by the following amount. 



(14) 



AacsB = a - a„„ = 1.6 ± 0.6 fm. 



(15) 
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Table 1. CSB differences of the ^So effective range parameters caused by nucleon 
mass splitting. 27r denotes the sum of all 27r-contributions and Trp the sum of all 
7r/9-contributions. TBE (non-iterative two-boson-exchange) is the sum of 2-k, np, 
and (ttct + ttuj). 





kin. en. 


OBE 


27r 




TTcr -t- iru! 


TBE 


Total 


Empirical 


^acsB (fm) 
^rcSB (fm) 


0.246 
0.004 


0.013 
0.001 


2.888 
0.055 


-1.537 
-0.031 


-0.034 
-0.001 


1.316 
0.023 


1.575 
0.028 


1.6 ±0.6 fm 
0.10 ±0.12 fm 



Recommended values for the corresponding effective ranges are [8] , 

r^„ = 2.75 ±0.11 fm, (16) 
r^p = 2.85 ± 0.04 fm, (17) 

implying 

ArcsB = <p-r^„ = 0.10 ± 0.12 fm. (18) 

Traditionally, it was believed that the meson mixing explains essentially all CSB 

effects. The largest contribution came from p — u mixing, and there was very meager 
knowledge of tt — r/ mixing. Recently, CSB was studied by the comparison of two 
charge-symmetric processes: D(7r+,r/)pp and 'D{'k~ ,rj)nn in the energy region of the 
rj threshold. The result for the ratio of the two processes in this energy region is 
R = da^ /da~ = 0.937±0.007, after a phase space correction is made for the difference 
in the threshold energies of the two reactions [21]. The deviation of R from 1 is an 
indication of CSB which is mostly due to tt — 77 mixing. A phenomenological fully 
relativistic model, which is based on coupled channel NTT-Nyy amplitudes, takes into 
account different nn and pp FSI and explicitly includes tt— 77 mixing, was developed [22] 
and compared to the data yielding for the tt — rj mixing angle the value of (1.5 ± 0.4) 
deg, consistent with the mixing angle determined from particle decays and isospin- 
forbidden processes as well as with several other theoretical predictions [21]. 

2.1.2. Theory. The difference between the masses of neutron and proton represents 
the most basic cause for CSB of the nuclear force. Therefore, it is important to have 
a very thorough accounting of this effect. 

The most trivial consequence of nucleon mass splitting is a difference in the kinetic 
energies: for the heavier neutrons, the kinetic energy is smaller than for protons. This 
raises the magnitude of the nn scattering length by 0.25 fm as compared to pp. 

Besides the above, nucleon mass splitting has an impact on all meson-exchange 
diagrams that contribute to the nuclear force. In 1998, the most comprehensive and 
thorough calculation of these CSB effects ever conducted has been published [23] . The 
investigation is based upon the Bonn Full Model for the NN interaction [24]. Here, we 
will summarize the results. For this we devide the total number of meson exchange 
diagrams that is involved in the nuclear force into several classes. Below, we report 
the results for each class. 

(i) One-boson-exchange (OBE, figure 1) contributions mediated by 7r'^(135), 
/9°(770), w(782), ao/(S(980), and cr'(550). In the Bonn Full Model [24], the 
a' describes only the correlated 2ti exchange in tttt — S'-wavc (and not the 
uncorrelated 27r exchange since the latter is calculated explicitly, cf. figure 2). 
Charge-symmetry is broken by the fact that for pp scattering the proton mass is 
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tt", a' ,p°, u}, 5 



P P 

(a) 



tt", a' ,p° ,u;, 5 



(b) 



Figure 1. One-boson-exchange (OBE) contributions to (a) pp and (b) nn 
scattering. 



used in the Dirac spinors representing the four external legs [figure f (a)], while 
for nn scattering the neutron mass is applied [figure 1(b)]. The CSB effect from 
the OBE diagrams is very small (cf. table 1). 

(ii) 27r-exchange diagrams. This class consists of three groups; namely the 
diagrams with NN, NA and AA intermediate states, where A refers to the baryon 
with spin and isospin | and mass 1232 MeV. The most important group is the 
one with A'' A intermediate states which wc show in figure 2. Part (a) of figure 2 
applies to pp scattering, while part (b) refers to nn scattering. When charged- 
pion exchange is involved, the intermediate-state nucleon differs from that of the 
external legs. This is one of the sources for CSB from this group of diagrams. 
The 2tt class of diagrams causes the largest CSB effect (cf. table 1 and dashed 
curve in figure 3). 

(iii) 7rp- exchanges. Graphically, the np diagrams can be obtained by replacing in 

each 2tt diagram (e. g., in figure 2) one pion by a p- meson of the same charge 
state. The effect is typically opposite to the one from 2n exchange. 

(iv) Further 37r and Aw contributions {wa + nLv). The Bonn potential also includes 
some 37r-exchanges that can be approximated in terms of ttct diagrams and 4ir- 
exchanges of ttlj type. The sum of the two groups is small, indicating convergence 
of the diagramatic expansion. The CSB effect from this class is essentially 
negligible. 

The total CSB difference of the singlet scattering length caused by nucleon mass 
splitting amounts to 1.58 fm (cf. table 1) which agrees well with the empirical value 
1.6 ± 0.6 fm. Thus, nucleon mass splitting alone can explain the entire empirical CSB 
of the singlet scattering length [25]. This is a remarkable result. 

The impact of the various classes of diagrams on CSB phase shift differences are 
shown in figure 3. The total effect is the largest in the ^So state where it is most 
noticable at low energy; e. g., at 1 MeV, the phase shift difference is 1.8 deg. The 
difference decreases with increasing energy and is about 0.15 deg at 300 MeV, in ^Sq. 

The CSB effect on the phase shifts of higher partial waves is small; in P and D 
waves, typically in the order of 0.1 deg at 300 MeV and less at lower energies. This 
fact may suggest that CSB in partial waves other than L = may be of no relevance. 
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Figure 2. Two-pion-cxchange contributions with ATA intermediate states to (a) 
pp and (b) nn scattering. 



In references [26] it was shown that this is not true: CSB beyond the S waves is crucial 
for the explanation of the Nolcn-Schiffcr anomaly. 

Before finishing this subsection, a word is in place concerning other mechanisms 
that cause CSB of the nuclear force. Traditionally, it was believed that —uj mixing 
explains essentially all CSB in the nuclear force [8]. However, recently some doubt 
has been cast on this paradigm. Some researchers [27, 28, 29, 30] found that p° — lo 
exchange may have a substantial dependence such as to cause this contribution to 
nearly vanish in NN. Our finding that the empirically known CSB in the nuclear force 
can be explained solely from nucleon mass splitting (leaving essentially no room for 
additional CSB contributions from — w mixing or other sources) fits well into this 
scenario. On the other hand, Miller [9] and Coon and coworkers [31] have advanced 
counter-arguments that would restore the traditional role of p-uj exchange. The issue 
is unresolved. Good summaries of the controversial points of view can be found in 
references [9, 32, 33]. 

Finally, for reasons of completeness, we mention that irreducible diagrams 
of TT and 7 exchange between two nucleons create a charge-dependent nuclear 
force. Recently, these contributions have been calculated to leading order in chiral 
perturbation theory [34]. It turns out that to this order the 717 force is charge- 
symmetric (but does break charge independence). 
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Figure 3. CSB phase shift differences 5nn — Spp (without electromagnetic 
interactions) for laboratory kinetic energies Ti^b below 300 MeV and partial waves 
with L < I. The CSB effects due to the kinetic energy, OBE, the entire 2n model, 
and TT/O exchanges are shown by the dotted, dash-triple-dot, dashed, and dash-dot 
curves, respectively. The solid curve is the sum of all CSB effects. 



2.2. Charge independence breaking 

The empirical values for the np singlet effective range parameters are [35]: 



a„p = -23.740 ± 0.020 fm, r„p = 2.77 ± 0.05 fm. (19) 
It is useful to define the following averages: 

a = ^{a^p + aZj = - 18.1 ± 0.6 fm, (20) 

r = Ik + r^n) = 2-80 ± 0.12 fm. (21) 
Ignoring CSB, the CIB differences in the effective range parameters are given by: 

AaciB =a-anp = 5.64 ± 0.60 fm, (22) 

Arc/B = f - r„p = 0.03 ± 0.13 fm. (23) 



The major cause of CIB in the NN interaction is pion mass splitting. Based upon 
the Bonn Full Model for the NN interaction [24], the CIB due to pion mass splitting 



Nucleon-nucleon interaction 



10 



p p 



n p p n 



P P 

(a) 



n p n p 

(b) 



Figure 4. One-pion exchange (OPE) contributions to (a) pp and (b) np 
scattering. 



has been calculated carefully and systematically in reference [36] . We will discuss now 
the various classes of diagrams and their contributions to CIB. 

(i) One-pion-exchange (OPE). The CIB effect is created by replacing the diagram 
figure 4(a) by the two diagrams figure 4(b). The effect caused by this replacement 
can be understood as follows. In nonrelativistic approximation| and disregarding 
isospin factors, OPE is given by 

V,A9.,m.) = S^^^^^^^^ F^r^Ni^^NN, |k|) (24) 

with M the average nucleon mass, ttItt the pion mass, and k the momentum 
transfer. The above expression includes a ttNN vertex form-factor, F-^nn, which 
depends on the cutoff mass A^jvtv and the magnitude of the momentum transfer 
|k|. For 5 = and T = 1, where S and T denote the total spin and isospin of the 
two-nucleon system, respectively, we have 

°^^i.(5.,m.) = ^^^^ ^ i^.V(A^iv;v, |k|) , (25) 

where the superscripts 01 refer to ST. In the state, this potential expression 
is repulsive. The charge-dependent OPE is then, 

oiyPf = 0Vi.(5.o,m.o) (26) 

t For pedagogical reasons, we use simple, aijproxiuiatc expressions to discuss the effects from 
pion exchange. Note, however, that in the calculations of reference [36] relativistic time-ordered 
perturbation theory is applied in its full complexity and without approximations. 



Table 2. CIB contributions to the ^Sq scattering length, AaciB, and effective 
range, ArciB, from various components of the NN interaction. 





OPE 


2tt 


7rp 


■KCT + TTLU 


Total 


Empirical 


^aciB (fm) 
^rciB (fm) 


3.243 
0.099 


0.360 
0.002 


-0.383 
-0.006 


1.426 
0.020 


4.646 
0.115 


5.64 ±0.60 
0.03 ±0.13 
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for pp scattering, and 

oiynp = 2 0Vi,(5.±,m,±) - 01^1,(5,0, m^o) (27) 

for np scattering. If we assume charge-independence of g„ (i. e., 5^0 = g.,^±), then 
all CIB comes from the charge splitting of the pion mass, which is [37] 

m^o = 134.977 MeV, (28) 
m,^± = 139.570 McV. (29) 

Since the pion mass appears in the denominator of OPE, the smaller 7r°-mass 
exchanged in pp scattering generates a larger (repulsive) potential in the ^Sq 
state as compared to np where also the heavier 7r^-mass is involved. Moreover, 
the Tr'^-exchange in np scattering carries a negative sign, which further weakens 
the np OPE potential. The bottom line is that the pp potential is more repulsive 
than the np potential. The quantitative effect on AaciB is such that it explains 
about 60% of the empirical value (cf. table 2). This has been know for a long 
time. 

Due to the small mass of the pion, OPE is a sizable contribution in all partial 
waves including higher partial waves; and due to the pion's relatively large mass 
splitting (3.4%), OPE creates relatively large charge-dependent effects in all 
partial waves (see dashed curve in figure 5). 

(ii) 27r-exchange diagrams. We now turn to the CIB created by the 27r exchange 
contribution to the NN interaction. There are many diagrams that contribute 
(see reference [36] for a complete overview). For our qualitative discussion here, 
we pick the largest of all 27r diagrams, namely, the box diagrams with NA 
intermediate states, figure 6. Disregarding isospin factors and using some drastic 
approximation, the amplitude for such a diagram is 

gt 72 f (i^p [<T-k5-k]2 



16M4 25 7 {2Tr)^ {ml + k-^yiEp + E^ -2Eg) 

X F^j,^{A^NN, |k|) F^^^iA^NA, |k|) , (30) 

where k = p — q with q the relative momentum in the initial and final state 
(for simplicity, we are considering a diagonal matrix element); Ep = 1/ + 
and E^ = y/ M\ + p2 with Ma = 1232 MeV the A-isobar mass; S is the spin 
transition operator between nucleon and A. For the ttNA coupling constant, 
/ttNA, the quark- model relationship f^^^ = H/Iatjv used [24]. 
For small momentum transfers k, this attractive contribution is roughly 
proportional to m~'^. Thus for the 2tt exchange, the heavier pions will provide 
less attraction than the lighter ones. Charged and neutral pion exchanges occur 
for pp as well as for np, and it is important to take the isospin factors carried by 
the various diagrams into account. They are given in figure 6 below each diagram. 
For pp scattering, the diagram with double tt^ exchange carries the largest factor, 
while double tt^ exchange carries only a small relative weight in np scattering. 
Consequently, pp scattering is less attractive than np scattering which leads to 
an increase of AaciB by 0.79 fm due to the diagrams of figure 6. The crossed 
diagrams of this type reduce this result and including all 27r exchange diagrams 
one finds a total effect of 0.36 fm [36]. 

(iii) 7rp-exchanges. This group is, in principle, as comprehensive as the 27r-exchanges 
discussed above. Graphically, the wp diagrams can be obtained by replacing in 
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Figure 5. CIB phase shift differences 5np — 5 [with S = {Spp + inn)/2] for 
laboratory Icinctic energies Ti^t below 300 MeV and partial waves with orbital 
angular momentum L < 1. The CIB effects due to OPE, the entire 2tt model, 
TTp exchanges, and (na + ww) contributions are shown by the dashed, dash-dot, 
dash-triple-dot, and dotted curves, respectively. The solid curve is the sum of all 
CIB effects. 



each 27r-diagram one of the two pions by a p-meson of the same charge-state. This 
contribution to CIB (dash-triple-dot curve in figure 5) is generahy small, and (in 
most states) opposite to the one from 2n. 

(iv) Further Stt and An contributions (ttct -f ttoj) . As discussed, the Bonn potential 
also includes some 37r-exchanges that can be approximated in terms of ira 
diagrams and 47r-exchanges of nu) type. These diagrams carry the same isospin 
factors as OPE. The CIB effect from this class is very small, except in ^So (dotted 
curve in figure 5). Notice that this effect has always the same sign as the effect 
from OPE (dashed curve), but it is substantially smaller. The reason for the OPE 
character of this contribution is that ttct prevails over nuj and, thus, determines 
the character of this contribution. Since sigma-exchange is negative and since, 
futhermore, the propagator in between the tt and the a exchange is also negative, 
the overall sign of the ttct exchange is the same as OPE. Thus, it is like a short- 
ranged OPE contribution. 
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Figure 6. 27r-exchange box diagrams with iVA intermediate states that 
contribute to (a) pp and (b) np scattering. The numbers below the diagrams 
are the isospin factors. 



Concerning the singlet scattering length, the CIB contributions discussed explain 
about 80% of AaciB (cf. table 2). Ericson and Miller [38] arrived at a very similar 
result using the meson-exchange model of Partovi and Lomon [39] . 

The sum of all CIB effects on phase shifts is shown by the solid curve in figure 5. 
Notice that the difference between the solid curve and the dashed curve (OPE) in 
that figure represents the sum of all effects beyond OPE. Thus, it is clearly seen that 
OPE dominates the CIB effect in all partial waves, even though there are substantial 
contributions besides OPE in some states, notably ^Sq and ^Pi. 

In reference [36], also the effect of rho-mass splitting on the ^Sq effective range 
parameters was investigated. Unfortunately, the evidence for rho-mass splitting is 
very uncertain, with the Particle Data Group [37] reporting m^o — mp± = 0.4 ± 0.8 
MeV. Consistent with this, nipO = 769 MeV and m^i = 768 MeV, i. e., a splitting 
of 1 MeV was assumed, in the exploratory study of reference [36]. With this, one 
finds AaciB = —0.29 fm from one-rho-exchange, and AaciB = 0.28 fm from the 
non-iterative np diagrams with NN intermediate states. Thus, individual effects are 
small and, in addition, there are substantial cancellations between the two classes 
of diagrams that contribute. The net result is a vanishing effect. Thus, even if the 
rho-mass splitting will be better known one day, it will never be a great source of CIB. 

Another CIB contribution to the nuclear force is irreducible pion-photon {ttj) 
exchange. Traditionally, it was believed that this contribution would take care of the 
remaining 20% of AaciB [38, 40, 41]. However, a recently derived wy potential based 
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upon chiral perturbation theory [34] decreases AaciB by about 0.5 fm, making the 
discrepancy even larger. 

Thus, it is a matter of fact that about 25% of the charge- dependence of the singlet 
scattering length is not explained — at this time. 

3. The ttA'^A'^ coupling constant 

For the nuclear force, the pion is the most important meson. Therefore, it is crucial 
to have an accurate understanding of the coupling of the pion to the nucleon. In the 
1990's, we have seen a controversial discussion about the precise value for the irNN 
coupling constant. Wc will first briefly review the events and then discuss in which 
way the NN data impose constraints on this important coupling constant. 

From 1973 to 1987, there was a consensus that the -kNN coupling constant is 
gl/A-K = 14.3 ± 0.2 (equivalent to fl = 0.079 ± 0.001§). This value was obtained 
by Bugg et al. [43] from the analysis of Tr^p data in 1973, and confirmed by Koch 
and Pietarinen [44] in 1980. Around that same time, the neutral-pion coupling 
constant was determined by KroU [45] from the analysis of pp data by means of 
forward dispersion relations; he obtained g^o/^T^ = 14.52 ± 0.40 (equivalent to 
/^o = 0.080 ±0.002). 

The picture changed in 1987, when the Nijmcgen group [46] determined the 
neutral-pion coupling constant in a partial-wave analysis of pp data and obtained 
5^0 /47r = 13.1 ± 0.1. Including also the magnetic moment interaction between protons 
in the analysis, the value shifted to 13.55 ± 0.13 in 1990 [47]. Triggered by these 
events, Arndt et al. [48] reanalysed the ir'^p data to determine the charged-pion 
coupling constant and obtained g^±/47r = 13.31 ± 0.27. In subsequent work, the 
Nijmegcn group also analysed np, pp, and ttN data [49]. The status of their work as 
of 1993 is summarized in Ref. [50] where they claim that the most accurate values 
are obtained in their combined pp and np analysis yielding g'^o/A-K = 13.47 ± 0.11 
(equivalent to fl„ = 0.0745 ± 0.0006) and gl±/ATT = 13.54 ± 0.05 (equivalent to 

= 0.0748 ± 0.0003). The latest analysis of all ix^p data below 2.1 GeV conducted 
by the VPI group using fixed-t and forward dispersion relation constraints has 
generated g^±/47r = 13.75 ± 0.15 [51]. The VPI NN analysis extracted g^o/47r « 13.3 
and (?^±/47r « 13.9 as well as the charge-independent value (?^/47r ~ 13.7 [52, 53]. 

Also Bugg and coworkers have performed new determinations of the -kNN 
coupling constant. Based upon precise ir^p data in the 100-310 MeV range and 
applying fixed-t dispersion relations, they obtained the value g'^^/in = 13.96 ± 0.25 
(equivalent to f^± = 0.0771 ± 0.0014) [54]. Prom the analysis of NN elastic data 
between 210 and 800 MeV, Bugg and Machleidt [55] have deduced gl±/4:n = 

§ Using ttNN Lagrangians as defined in the authoritative review [42], the relevant relationships 
between the pseudoscalar pion coupling constant, g^, and the pseudovector one, /x, are 



2 



47r V ni^± J 

and 



f-^o = 180.773/^0 (31) 



9l± f Mp + Mn'^^ 



47r V rn^± 



181.022/;^± . (32) 



with Mp = 938.272 MeV the proton mass, M„ = 939.566 MeV the neutron mass, and m^± = 139.570 
MeV the mass of the charged pion. 
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Table 3. Important coupling constants and the predictions for the deuteron and 
some pp phase shifts for five models discussed in the text. 





A 


L! C D 


E 


Eiii|)i].ical 




Important coupling constants 








13.6 


13.6 14.0 14.4 


13.6 






13.6 


13.6 14.0 14.4 


14.4 






6.1 


3.7 6.1 6.1 


6.1 








The deuteron 






Q (fm2) 


0.270 


0.278 0.276 0.282 


0.278 


0.276(2)" 


V 


0.0255 


0.0261 0.0262 0.0268 


0.0264 


0.0256(4)'' 


As (fm-l/2) 


0.8845 


0.8842 0.8845 0.8845 


0.8847 


0.8845(8)= 


Pd (%) 


4.83 


5.60 5.11 5.38 


5.20 








^Po pp phase shifts (dcg) 






10 MeV 


3.726 


4.050 3.881 4.039 


3.726 


3.729(17)'' 


25 MeV 


8.588 


9.774 8.981 9.384 


8.588 


8.575(53)'' 


50 MeV 


11.564 


14.070 12.158 12.763 


11.564 


11.47(9)=* 



" Corrected for meson-exchange currents and relativity. 
Reference [62]. 
Reference [63]. 

Nijmegen pp multi-energy phase shift analysis [64]. 



13.69 ± 0.39 and .g2o/47r = 13.94 ± 0.24. 

Thus, it may appear that recent determinations show a consistent trend towards 
a lower vahie for g^r with no indication for substantial charge dependence. 

However, tliis is not true and for a comprehensive overview of recent 
determinations of the ttNN coupling constant, see reference [56]. In particular, 
there is one determination that does not follow the above trend. Using a modified 
Chew extrapolation procedure, the Uppsala Neutron Research Group has deduced 
the charged-pion coupling constant from high precision np charge-exchange data at 
162 MeV [57]. Their latest result is 5^±/47r = 14.52 ± 0.26 [58]. We note that the 
method used by the Uppsala Group is controversial [59, 60]. 

Since the pion plays a crucial role in the creation of the nuclear force, many NN 
observables are sensitive to the 7rA''A'' coupling constant, Qt^. We will discuss here the 
most prominent cases and their implications for an accurate value of ^tt- 

We will focus on the deuteron, NN analyzing powers Ay, and the singlet scattering 
length. Other NN observables with sensitivity to g^^ are spin transfer coefficients. 
Concerning the latter and their implications for g^^, we refer the interested reader to 
references [55, 61]. 

3.1. The deuteron 

The crucial deuteron observables to consider are the quadrupole moment, Q, and the 
asymptotic D/S state ratio, r/. The sensitivity of both quantities to g,^ is demonstrated 
in table 3. The calculations are based upon the CD-Bonn potential [65, 66]) which 
belongs to the new generation of high-precision NN potentials that fit the NN data 
below 350 MeV with a 'perfect' x^/datum of about one. The numbers in table 3 are 
an update of earlier calculations of this kind [67, 68] in which older NN potentials 
were applied. However, there are no substantial differences in the results as compared 
to the earlier investigations. 
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Figure 7. ^Po phase shifts of proton-proton scattering as predicted by Model A 
and E (g^o/^T = 13.6, solid line), B (k^ = 3.7, dash-3dot), C (ff^o/^Tr = 14.0, 
dashed), and D (g^g/47r = 14.4, dash-dot) ||. The solid dots represent the 
Nijmegen pp multi-energy phase shift analysis [64]. 



For meaningful predictions, it is important that all deuteron models considered are 

realistic. This requires that besides the deuteron binding energy (that is accurately 
reproduced by all models of table 3) also other empirically well-known quantities 
are correctly predicted, like the deuteron radius, r^, and the triplet effective range 
parameters, at and rt. As it turns out, the latter quantities are closely related to the 
asymptotic S-state of the deuteron, As, which itself is not an observable. However, it 
has been shown [63] that for realistic values of rj,, at, and rt, the asymptotic S-state 
of the deuteron comes out to be in the range As = 0.8845 ± 0.0008 fm^i/^ xhus. 
As plays the role of an important control number that tells us if a deuteron model is 
realistic or not. As can be seen from table 3, all our models pass the test. 

Model A of table 3 uses the currently fashionable value for the ttNN coupling 
constant g'^/'in = 13.6 which clearly underpredicts Q while rj is predicted satisfactorily. 
One could now try to fix the problem with Q by using a weaker p-meson tensor- 
coupling to the nucleon, fp. It is customary to state the strength of this coupling in 
terms of the tensor-to- vector ratio of the p coupling constants, Kp = fp/gp- Model 
A uses the 'large' value Kp = 6.1 recommended by Hoehler and Pietarinen [69]. 
Alternatively, one may try the value implied by the vector-meson dominance model 
for the electromagnetic form factor of the nucleon [70] which is Kp — 3.7. This is done 
in our Model B which shows the desired improvement of Q. However, a realistic model 
for the NN interaction must not only describe the deuteron but also NN scattering. 
As discussed in detail in reference [71], the small Kp cannot reproduce the ei mixing 
parameter correctly and, in addition, there are serious problems with the ''Pj phase 
shifts, particularly, the ^Pq (cf. lower part of table 3 and figure 7). Therefore, Model 
B is unrealistic and must be discarded. 

The only parameters left to improve Q are and the t:NN vertex form-factor, 
FnNN (cf. equation 24, above). As for the p meson, F^nn is heavily constrained by 
NN phase parameters, particularly, ei. The accurate reproduction of ei as determined 
in the Nijmegen np multi-energy phase shift analysis [64] essentially leaves no room 
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Table 4. x'^/datum for the fit of the world Ay data below 350 MeV (subdivided 
into three energy ranges) using different values of the nNN coupling constant. 





Coupling constant 


9^0 /47r 


Energy range (# of data) 


13.2 


13.6 


14.0 


14.4 






A 


C 


D 


0-17 MeV (45 data) 


0.84 


1.43 


2.71 


4.66 


17-125 MeV (148 data) 


1.05 


1.06 


1.54 


2.45 


125-350 MeV (624 data) 


1.24 


1.22 


1.26 


1.34 



for variations of Fj^nn once the p meson parameters are fixed. 

Thus, we are finally left with only one parameter to fix the Q problem, namely 
Qn. As it turns out, for relatively small changes of Q^/^tt there is a linear relationship, 
as demonstrated in table 3 by the predictions of Model A, C and D which use 
g'^/^TT = 13.6, 14.0, and 14.4, respectively. Consistent with earlier studies [67, 68], one 
finds that 5^/47r > 14.0 is needed to correctly reproduce Q. 

However, a pion coupling with g^/^Tr > 14.0 creates problems for the ^Pq phase 
shifts which are predicted too large at low energy (cf. lower part of table 3 and figure 7). 
Now, a one-boson-exchange (OBE) model for the NN interaction includes several 
parameters (about one dozen in total). One may therefore try to improve the ^Pq by 
readjusting some of the other model parameters. The vector mesons (p and lo) have 
a strong impact on the ^Pq (and the other P waves). However, due to their heavy 
masses, they are more effective at high energies than at low ones. Therefore, p and w 
may produce large changes of the ^Pq phase shifts in the range 200-300 MeV, with little 
improvement at low energies. The bottom line is that in spite of the large number of 
parameters in the model, there is no way to fix the ^Pq phase shift at low energies. In 
this particular partial wave, the pion coupling constant is the only effective parameter, 
at energies below 100 MeV. The pp phase shifts of the Nijmegen analysis [64] as well 
as the pp phases produced by the VPI group [72] require g'^/An < 13.6. 

Notice that this finding is in clear contradiction to our conclusion from the 
deuteron Q. 

There appears to be a way to resolve this problem. One may assume that 
the neutral pion, tt*^, couples to the nuclcon with a slightly different strength than 
the charged pions, tt^. This assumption of a charge-splitting of the ttA^A^ coupling 
constant is made in our Model E where WG USG Q^o 

/An = 13.6 and gl±/4:Tr = 14.4. 
This combination reproduces the pp ^Pq phase shifts at low energy well and creates a 
sufficiently large deuteron Q. 

3.2. Analyzing powers 

In our above considerations, some pp phase shifts played an important role. In 
principle, phase shifts are nothing else but an alternative representation of data. Thus, 
one may as well use the data directly. Since the days of Gammel and Thaler [4], it 
is well-known that the triplet P-wave phase shifts arc fixed essentially by the NN 
analyzing powers, Ay. Therefore, we will now take a look at Ay data and compare 
them directly with model predictions. 

In figure 8, we show high-precision pp Ay data at 9.85 MeV from Wisconsin [73]. 
The theoretical curves shown are obtained with g^o/An = 13.2 (dotted), 13.6 (solid). 
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Figure 8. The proton-proton analyzing power Ay at 9.85 MeV. Tiie tlieoretical 
curves are calculated with g^o/iir = 13.2 (dotted), 13.6 (solid, Model A), and 
14.4 (dash-dot, Model D) and fit the data with a x^/datum of 0.98, 2.02, and 
9.05, respectively. The solid dots represent the data taJcen at Wisconsin [73]. 



and 14.4 (dash-dot) and fit the data with a x^/datum of 0.98, 2.02, and 9.05, 
respectively. Clearly, a small coupling constant around 13.2 is favored. Since a single 
data set is not a firm basis, we have looked into all pp Ay data in the energy range 
0-350 MeV. Our results are presented in table 4 where we give the /datum for the 
fit of the world pp Ay data below 350 MeV (subdivided into three energy ranges) for 
various choices of the neutral -kNN coupling constant. It is seen that the pp Ay data 
at low energy, particularly in the energy range 0-17 MeV, are very sensitive to the 
ttA^'A'' coupling constant. A value g^o/^i^ < 13.6 is clearly preferred, consistent with 
what we extracted from the single data set at 9.85 MeV as well as from the ^Pq phase 
shifts. 

Next, we look into the np Ay data. A single sample is shown in figure 9, the np 
Ay data at 12 MeV from TUNL [74]. Predictions are shown for Model A (sohd line), 
D (dash-dot), and E (dash-triplc-dot). The charge-splitting Model E fits the data best 
with a x^/datum of 1.00 (cf. table 5). We have also considered the entire np Ay data 
measured by the TUNL group [74] in the energy range 7.6-18.5 MeV (31 data) as well 
as the world np Ay data in the energy ranges 17 MeV (120 data). It is seen that 
there is some sensitivity to the TriVA'' coupling constant in this energy range, while 
there is little sensitivity at energies above 17 MeV (cf. table 5). 

Consistent with the trend seen in the 12 MeV data, the larger data sets below 17 
MeV show a clear preference for a coupling constant around 14.4 if there is no charge 
splitting of g^Tr- This implies that without charge-splitting it is impossible to obtain 
an optimal fit of the pp and np Ay data. To achieve this best fit, charge-splitting is 
needed, like (?^o/47r = 13.6 and g^±/A-K = 14.0, as considered in column 5 of table 5. 
The drastic charge-splitting of Model E is not favored by the more comprehensive np 
Ay data sets. 

The balance of the analysis of the pp and np Ay data then is: 5^0/477 < 13.6 
and g^^/A-K > 14.0. Notice that this splitting is consistent with our conclusions from 
the deuteron. Thus, we have now some indications for charge-splitting of from two 
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Figure 9. The neutron-proton analyzing power Ay at 12 McV. The tlicorctical 
curves arc calculated with g^Q/4n = g'^^/4n = 13.6 (solid line, Model A), 
9^0 Z^"" = = 14.4 (dash-dot, Model D), and the charge-splitting p^(,/47r = 

13.6, g^^/47r = 14.4 (dash-3dot. Model E). The soUd dots represent the data 
taken at TUNL [74]. 

very different observables, namely the deuteron quadrupole moment and np analyzing 
powers. 

Therefore, it is worthwhile to look deeper into the issue of charge-splitting of 
the ttNN coupling constant. Unfortunately, there are severe problems with any 
substantial charge-splitting for two reasons. First, theoretical work [77] on isospin 
symmetry breaking of the ttNN coupling constant based upon QCD sum rules comes 
up with a splitting of less than 0.5% for and, thus, cannot explain the large charge 
splitting indicated above. Second, a problem occurs with the conventional explanation 
of the charge-dependence of the singlet scattering length, which we will explain now. 



Table 5. x'^ /datum for the fit of various sets of np Ay data using different values 
for the ttNN coupling constants. 

Coupling constants g^g/iTr; 9^±/47r 



Energy, data set (# of data) 


13.6; 13.6 


14.0; 14.0 


14.4; 14.4 


13.6; 14.0 


13.6; 14.4 




A 


C 


D 




E 


12 MeV [74] (9 data) 


2.81 


2.27 


1.79 


1.53 


1.00 


7.6-18.5 MeV [74] (31 data) 


1.89 


1.56 


1.29 


1.28 


1.32 


0-17 MeV world data (120) 


1.17 


1.03 


0.94 


0.99 


1.19 


17-50 MeV [75] (85 data) 


1.16 


1.12 


1.14 


1.18 


1.18 


17-125 MeV world data (416) 


0.89 


0.89 


0.91 


0.91 


0.94 
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3.3. Charge- dependence of the singlet scattering length and charge-dependence of the 
pion coupling constant 

Here, we are going to show in detail how charge-sphtting of the ttNN couphng constant 
affects the charge-dependence of the scattering length. It will turn out that the 
suggested charge-splitting of g„ causes a disaster for our established understanding of 
the charge-dependence of the singlet scattering length. 

Our above considerations suggest charge-splitting of g,r, like 

glo/^TT = 13.6 , (33) 
5^±/47r=14.4, (34) 
cf. Model E of table 3. We will now discuss how this charge-splitting of gi^ affects 
Aac/B (more details can be found in the original paper reference [76]). 

Accidentally, this splitting is — in relative terms — about the same as the pion-mass 
splitting; that is 

« ^ . (35) 

As discussed (cf. equations (25) and (30) and text below these equations), for zero 
momentum transfer, we have roughly for one-pion exachange 

OPE^f^^ (36) 

\ rn^ J 

and for 27r exchange 

TPE ~ ( ^] , (37) 

which is not unexpected, anyhow. On the level of this qualitative discussion, we can 
then predict that any pionic charge-splitting satisfying equation (35) will create no GIB 
from pion exchanges. Consequently, a charge-splitting of g.^ as given in equations (33) 
and (34) will wipe out our established explanation of GIB of the NN interaction. 

In reference [76], accurate numerical calculations based upon the Bonn meson- 
exchange model for the NN interaction [24] have been conducted. The details of these 
calculations are spelled out in reference [36] where, however, no charge-splitting of g^^ 
was considered. Assuming the of equations (33) and (34), one obtains the AaciB 
predictions given in the last column of table 6. It is seen that the results of an accurate 
calculation go even beyond what the qualitative estimate suggested: the conventional 
GIB prediction is not only reduced, it is reversed. This is easily understood if one 
recalls [cf. equations (25) and (30)] that the pion mass appears in the propagator 
(m^ -|- k^)~^. Assuming an average « m^, the 7% charge splitting of will 
lead to only about a 3% charge-dependent effect from the propagator. Thus, if a 6% 
charge-splitting of g^ is used, this will not only override the pion-mass effect, it will 
reverse it. 

Based upon this argument and on the numerical results, one can then estimate 
that a charge-splitting of 5 J of only about 3% (e. g-,5^o/47r= 13.6 and 5^±/47r = 14.0) 
would erase all predictions of GIB in the singlet scattering length derived from pion 
mass splitting. 

Besides pion mass splitting, we do not know of any other essential mechanism 
to explain the charge-dependence of the singlet scattering length. Therefore, it is 
unlikely that this mechanism is annihilated by a charge-splitting of g^r . This may be 
taken as an indication that there is no significant charge splitting of the TriVA^ coupling 
constant. 
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Table 6. Predictions for AaciB in units of fm without and with the assumption 
of charge-dependence of Qtt ■ 





No charge-dependence of 


Charge-dependent g-,r' 




9^o/47r = 92^/47r = 14.4 


/4n = 13.6 






(,2^/47r = 14.4 


iTT 


3.24 


-1.58 


2n 


0.36 


-1.94 


wp, wcr, 1TU> 


1.04 


-0.97 


Sum 




l.il) 


Empirical 


5.64 ± 0.60 



3.4- Conclusions 

Several NN observables can be identified that are very sensitive to the wNN coupling 
constant, g^^. They all carry the potential to determine g^r with high precision. 

In particular, we have shown that the pp Ay data below 17 MeV are very sensitive 
to and imply a value g^/Aw » 13.2. The np Ay data below 17 MeV show moderate 
sensitivity and the deuteron quadrupole moment shows great sensitivity to gj^; both 
np observables imply g'^/Air > 14.0. 

The two different values may suggest a relatively large charge-splitting of g^. 
However, a charge-splitting of this kind would completely destroy our established 
explanation of the charge-dependence of the singlet scattering length. Since this is 
unlikely to be true, we must discard the possibility of any substantial charge-splitting 
of 5^. 

The conclusion then is that we are faced with real and substantial discrepancies 
between the values for g,^ based upon different NN observables. The reason for this 
can only be that there are large, unknown systematic errors in the data and/or large 
uncertainties in the theoretical methods. Our homework for the future is to find these 
errors and eliminate them. 

Another way to summarize the current confused situation is to state that, 
presently, any value between 13.2 and 14.4 is possible for (7^/47r depending on which 
NN observable you pick. If we want to pin down the value more tightly, then we are 
faced with three possible scenarios: 

• g^ is small, g^/4:TT < 13.6: 

The deuteron rj and pp scattering at low energies are described well; there 
are moderate problems with the np Ay data below 17 MeV. The most serious 
problem is the deuteron Q. Meson-exchange current contributions (MEG) and 
relativistic corrections for Q of 0.016 fm^ or more would solve the problem. 
Present calculations predict about 0.010 fm^ or less. A serious reinvestigation 
of this issue is called for. 

• 9n is large, gl/^n > 14.0: 

The deuteron Q is well reproduced, but r] is predicted too large as compared 
to the most recent measurement by Rodning and Knutsen [62], r] = 0.0256(4). 
Note, however, that all earlier measurements of rj came up with a larger value; 
for example, Borbely et al. [78] obtained t] = 0.0273(5). There are no objectively 
verifiable reasons why the latter value should be less reliable than the former one. 
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The deuteron rj carries the potential of being the best observable to determine g^^ 

(as pointed out repeatedly by Ericson [63, 79] in the 1980's); but the unsettled 
experimental situation spoils it all. The np Ay data at low energy are described 
well. The most serious problem are the pp Ay data below 100 MeV. 

• is 'in the middle', 13.6 < gl/An < 14.0: 

we have all of the above problems, but in moderate form. 

In conclusion, to arrive at an accurate value for g^, there is a lot of homework to 
do — for theory and experiment. 

4. Phase shift analysis 

In spite of the large NN database available in the 1990's, conventional phase shift 
analyses are by no means perfect. For example, the phase shift solutions obtained by 
Bugg [80] or the VPI/GWU group [72] typically have a x^/datum of 1.3 or more, for 
the energy range 0-425 MeV. This may be due to inconsistencies in the data as well 
as deficiencies in the constraints applied in the analysis. In any case, it is a matter 
of fact that within the conventional phase shifts analysis, in which the lower partial 
waves are essentially unconstrained, a better fit cannot be achieved. 

About two decades ago, the Nijmegen group embarked on a program to 
substantially improve NN phase shift analysis. To achieve their goal, the Nijmegen 
group took two decisive measures [64]. First, they 'pruned' the database; i.e., they 
scanned very critically the world NN database (all data in the energy range 0-350 MeV 
laboratory energy published in a regular physics journal between January 1955 and 
December 1992) and eliminated all data that had either an improbably high (more 
than three standard deviations off) or an improbably low x^; of the 2078 world pp 
data below 350 MeV 1787 survived the scan, and of the 3446 np data 2514 survived. 
Second, they introduced sophisticated, semi-phenomenological model assumptions into 
the analysis. Namely, for each of the lower partial waves ( J < 4) a different energy- 
dependent potential is adjusted to constrain the energy-dependent analysis. Phase 
shifts are obtained using these potentials in a Schroedinger equation. From these 
phase shifts the predictions for the observables are calculated including the x^ for the 
fit of the experimental data. This x^ is then minimized as a function of the parameters 
of the partial-wave potentials. Thus, strictly speaking, the Nijmegen analysis is a 
potential analysis; the final phase shifts are the ones predicted by the 'optimized' 
partial-wave potentials. 

In the Nijmegen analysis, each partial-wave potential consists of a short- and a 
long-range part, with the separation line at r — 1.4 fm. The long-range potential Vl 
(r > 1.4 fm) is made up of an electromagnetic part Vem and a nuclear part Vat: 



for neutron-proton scattering, where Vc denotes an improved Coulomb potential 
(which takes into account the lowest-order relativistic corrections to the static 



Vl = Vem + Vn 



(38) 



The electromagnetic interaction can be written as 
Vem{pp) = Vc + Vvp + Vmm{pp) 

for proton-proton scattering and 

VEM{np) = VMM{np) 



(40) 



(39) 
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Coulomb potential and includes contributions of all two-photon exchange diagrams); 

Vvp is the vacuum polarization potential, and Vmm the magnetic moment interaction. 

The nuclear long-range potential Vn consists of the local one-pion-exchange 
(OPE) tail ViTT (the coupling constant g-,^ being one of the parameters used to minimize 
the x^) multiplied by a factor M/E and the tail of the hcavy-boson-exchange (HBE) 
contributions of the Nijmegen78 potential [81] Vhbe, enhanced by a factor of 1.8 in 
singlet states; i. e. 

M 

V^AT = X Vi^{g^,m^) + f{S) x Vhbe (41) 

with f{S = 0) = 1.8 and f{S = 1) = 1.0, where S denotes the total spin of the 
two-nucleon system. The energy-dependent factor M/E (with E = ^^/A/P + q^, = 
AITiab/2) takes into account relativity in a 'minimal' way, damping the nonrelativistic 
OPE potential at higher energies. 

As indicated, Vi„ depends on the ttA^'A'' coupling constant and the pion mass 
m„, which gives rise to charge dependence. For pp scattering, the OPE potential is 

Vr: = V^^{g^o,m„o) (42) 

with nij^o the mass of the neutral pion. In np scattering, we have to distinguish 
between T = 1 and T = 0: 

Vr/iT) = -Fi.(5.o,m,o) + {-lf+^2Vi^{g^±,m^±) (43) 

The partial- wave short-range potentials (r < 1.4 fm) are energy-dependent 
square- wells (see figures 2 and 3 of reference [64]). The energy-dependence of the 
depth of the square- well is parametrized in terms of up to three parameters per partial 
wave. For the states with J < 4, there are a total of 39 such parameters (21 for pp 
and 18 for np) plus the pion-nucleon coupling constants (g^o and g^±). 

In the Nijmegen np analysis, the T = 1 np phase shifts are calculated from 
the corresponding pp phase shifts (except in ^Sq where an independent analysis 
is conducted) by applying corrections due to electromagnetic effects and charge 
dependence of OPE. Thus, the np analysis determines ^So{np) and the T = states, 
only. 

In the combined Nijmegen pp and np analysis [64], the fit for 1787 pp data and 
2514 np data below 350 MeV, available in 1993, results in the 'perfect' x^/datum = 
0.99. 

5. Recent additions to the NN database 

The world NN data (below 350 MeV) published before December 1992 have been listed 
and analyzed carefully by the Nijmegen group in their papers about the Nijmegen 
phase shift analysis [47, 64]. Therefore, we will focus here on data published after 
1992. 

5.1. Proton-proton data 

In the past decade, there has been a major breakthrough in the development of 
experimental methods for conducting hadron-hadron scattering experiments. In 
particular, the method of internal polarized gas targets applied in stored, cooled beams 
is now working perfectly in several hadron facilities, e. g., lUCF (Indiana, USA) and 
COSY (Jiilich, Germany). Using this new technology, lUCF has produced a large 
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Table 7. After-1992 proton-proton data below 350 MeV. 'Error' refers to the 
normalization error. This table contains 1127 observables and 32 normalizations 
resulting in a total of 1159 data. 



Tiab (McV) # observable Error (%) Institution(s) Ref. 



SnO-O 407 


14 a 


None 


Miinstcr 


f82l 


25.68 


8 D 


1.3 


Erlangen 


, Ziirich, PSl 


f83l 


25.68 


6 R 


1.3 


Erlangcn 


, Ziirich, PSl 


f83l 


25.68 


2 A 


1.3 


Erlangen 


, Ziirich, PSl 


f83l 


197.4 


41 P 


1.3 


Wisconsin, lUCF 


f84l 


197.4 


41 4t.t. 


2.5 


Wisconsin, lUCF 


f84l 


197.4 


41 A-,,-,, 


2.5 


Wisconsin, lUCF 


[841 


197.4 


41 Axz 


2.5 


Wisconsin, lUCF 


f841 


197.4 


39 Azz 


2.0 


Wisconsin, lUCF 


f85l 


197.8 


14 P 


1.3 


Wisconsin, lUCF 


f86l 


197.8 


14 A--,. 


2.4 


Wisconsin, lUCF 


f86l 

[OUJ 


197.8 


14 A.,., 


2.4 


Wisconsin, lUCF 


f86l 

[OUJ 


197.8 


14 A^^ 

^XZ 


2.4 


Wisconsin, lUCF 


[OUJ 


197.8 


10 D 




lUCF 


fell 


197.8 


5 R 


None 


lUCF 


fell 

[U -LJ 


197.8 


5 R' 


None 


lUCF 


fell 


197.8 


5 A 


None 


lUCF 


fell 

[U±J 


197.8 


5 A' 


None 


lUCF 


fell 

[U -LJ 


250.0 


41 P 


1.3 


lUCF, 


Wisconsin 


f87l 

[O 1 J 


250.0 


41 /Itt 

^ J- ^^xx 


2.5 


lUCF, 


Wisconsin 


f87l 

[O i J 


250.0 


41 A.,., 


2.5 


lUCF, 


Wisconsin 


f87l 

[O i J 


250.0 


41 Axz 


2.5 


lUCF, 


Wisconsin 


[87] 


280.0 


41 P 


1.3 


lUCF, 


Wisconsin 


[87] 


280.0 


41 Axx 


2.5 


lUCF, 


Wisconsin 


[87] 


280.0 


41 Ayy 


2.5 


lUCF, 


Wisconsin 


[87] 


280.0 


41 A:,^ 


2.5 


lUCF, 


Wisconsin 


[87] 


294.4 


40 P 


1.3 


lUCF, 


Wisconsin 


[87] 


294.4 


40 A^^ 


2.5 


lUCF, 


Wisconsin 


[87] 


294.4 


40 Ayy 


2.5 


lUCF, 


Wisconsin 


[87] 


294.4 


40 


2.5 


lUCF, 


Wisconsin 


[87] 


310.0 


40 P 


1.3 


lUCF, 


Wisconsin 


[87] 


310.0 


40 Axx 


2.5 


lUCF, 


Wisconsin 


[87] 


310.0 


40 Ayj, 


2.5 


lUCF, 


Wisconsin 


[87] 


310.0 


40 Axz 


2.5 


lUCF, 


Wisconsin 


[87] 


350.0 


40 P 


1.3 


lUCF, 


Wisconsin 


[87] 


350.0 


40 Aa^a; 


2.5 


lUCF, 


Wisconsin 


[87] 


350.0 


40 Ayy 


2.5 


lUCF, 


Wisconsin 


[87] 


350.0 


40 


2.5 


lUCF, 


Wisconsin 


[87] 



number of pp spin correlation parameters of very high precision. In table 7, we list 
the new lUCF data together with other pp data below 350 MeV published between 
January 1993 and December 1999. The total number of after-1992 pp data is 1159, 
which should be compared to the number of pp data in the (Nijmegen) 1992 base, 
namely, 1787. Thus, the pp database has increased by about 2/3 since 1992. The 
importance of the new pp data is further enhanced by the fact that they are of much 
higher quality than the old ones. 

The /datum produced by some recent phase shift analyses (PSA) and NN 
potentials in regard to the old and new databases are given in table 8. In this table, 
the '1992 database' is the Nijmegen database [47, 64] and the '1999 database' is the 
sum of the 1992 base and the after-1992 data. 
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Table 8. x^/datum for the NN data below 350 MeV applying some recent phase 
shift analyses (PSA) and NN potentials. 





VPl/GWU 
PSA [88] 


Nijmegen 1993 
PSA [64] 


Argonne Vig 
pot. [89] 


CD-Bonn 
pot. [66] 




proton-proton data 






1992 pp database (1787 data) 


1.28 


1.00 


1.10 


1.00 


After-1992 pp data (1145 data)" 


1.08 


1.24 


1.74 


1.03 


1999 pp database (2932 data)" 


1.21 


1.09 


1.35 


1.01 




neutron-proti 


on data 






1992 np database (2514 data) 


1.19 


0.99 


1.08 


1.03 


After-1992 np data (544 data)*" 


o.gs'^ 


0.99 


1.02 


0.99 


1999 np database (3058 data)*" 


1.16^= 


0.99 


1.07 


1.02 



"Without the 14 pp a data of Rcf. [82]. 
'Without after- 1992 np cr data, except [92]. 
■^Without the data of reference [91]. 



What stands out in table 8 are the rather large values for the /datum generated 
by the Nijmegen analysis and the Argonne potential for the the after-1992 pp data, 
which are essentially the new lUCF data. This fact is a clear indication that these new 
data provide a very critical test/constraint for any NN model. It further indicates that 
fitting the pre- 1993 pp data does not nessarily imply a good fit of those lUCF data. On 
the other hand, fitting the new lUCF data does imply a good fit of the pre-1993 data. 
The conclusion from these two facts is that the new lUCF data provide information 
that was not contained in the old database. Or, in other words, the pre-1993 data 
were insufficient and still left too much lattitude for pinning down NN models. One 
thing in particular that we noticed is that the phase shifts above 100 MeV have 
to be lower than the values given in the Nijmegen analysis. 

5.2. Neutron-proton data 

Neutron-proton data published between January 1993 and December 1999 are listed 
in table 9. Particular attention deserve the TUNL data on Actl and Actt between 5 
and 20 MeV [91] which made it possible to pin down the ei mixing parameter with 
unprecedented precision. 

Table 9 includes several new measurements of np differential cross sections (cr), 
with the largest sets produced by the Freiburg group [100] (871 data between 200 and 
350 MeV), the Uppsala group [105, 98] (92 and 162 MeV, 109 data), and at Louvain- 
la-Neuve [96] (84 data between 29 and 73 MeV). In table 10 we show the /datum 
for the reproduction of these data sets by some recent PSA and NN potentials. One 
observes that none of the PSA and potentials can reproduce these data accurately. 
For comparison, table 10 includes the x^/datum for the np a data by Bonner and 
coworkers [104] (652 data between 162 and 344 MeV, published in 1978) which are 
well reproduced. The large differences in the x^ implies that there are inconsistencies 
in the data. Since the PSA and potentials of table 10 were fitted to a database that 
includes the Bonner data but excludes the Freiburg data, one might suspect that 
this could be the explanation of the good x^ for the Bonner data and the bad one for 
Freiburg. Arndt has investigated this question [106] and found that this is not entirely 
true. When he excludes the Bonner data from the VPI/GWU analysis and uses the 
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Table 9. Aftcr-1992 iicutron-yjrotoii data below 350 McV. 'Error' refers to the 
noriiiaiizalioii error. 



Tiab (MeV) 


# observable 


Error (%) 


Institution (s) 


Ref. 


Q 1 1 K 


Q A /-T--^ 

y l\(tt 


None 




ronl 
[yuj 


4.ycs— iy. ( 




None 




roi 1 


A QSJ 1 ? 1 


A ^ 


None 




Fqi 1 
[91 J 


1 y1 1 1 
14. li 


O (T 


U. / 


Tiibingen 


Fool 

[yzj 


1 cr Q 
io.o 




None 


Bonn 


FriQl 
[yoj 




1 A ,-1- 


None 


Prague 


[y4J 


1 O 


1 A ^ 


None 


Prague 


Fnc;l 

[yoj 


9Q n 


yj u 


A n 


Lou vain- la-Neuve 




oi.o 


D cr 


4.U 


Louvain-la-Neuve 


FOftl 

[yoj 


O'4:.0 


O CT 


A n 


Lou vain- la-Neuve 


[yoj 


O / .0 


ft ^ 
O CT 


/I n 
4.U 


Louvain-la-Neuve 


FOftl 

[yoj 


4i.U 


CT 


/I n 

4.U 


Louvain-la-Neuve 


FOftl 

[yoj 


/I n 

40. U 


ft ^ 
O CT 


/I n 

4.U 


Louvain-la-Neuve 


FOftl 

[yoj 


/in n 
4y.U 


ft ^ 
D CT 


4.U 


L ou vain- 1 a- Neuve 


Fnftl 

[yoj 


Oo.U 


D CT 


A n 


Louvain-la-Neuve 


L'^'-'l 


Oo.U 


ft rt- 
D CT 


/t n 

4.U 


Lou vain- la-Neuve 


FQftl 

[yoj 


ox. 1 


ft rt- 


A n 

4.U 


Lou vain- la-Neuve 


[yoj 


D 1 . ( 


i-O CT 


r ioat 


Dasei, x^oi 


rQ7i 
[ycj 


fl7 7 


ft rt- 
D CT 


/t n 

4.U 


Louvain-la-Neuve 


FQftl 

[yoj 


70 8 


ft rt- 


A n 

4.U 


Louvain-la-Neuve 


FOftl 

[yoj 




04 CT 


z.o 


Uppsala 


FqqI 

[yoj 


1 7c^ OK 


o4 J 


4.y 


1 rtiUIVlr 


Fool 

[yyj 


1 QQ Q 


1 n9 /T 


^ n 


rreiDurg, roi 


Fl nnl 

[lUUJ 


ZUO. lO 


1 nn p 


A 7 


TPTTTA/TP 


[QQl 

[yyj 


917 OA 


1 nn p 


^.0 


1 ixi u ivrr 


[yyj 


ziy.o 


1U4 (J 


o.U 


Freiburg, PSI 


F l nnl 
[lUUJ 


0/in o 

Z4U.Z 


1 n7 >-r 


o.U 


Freiburg, PSI 


\^ nnl 

[lUUJ 


zDU.U 


Si R. 


o.U 




Ft m 1 
[lUlJ 


zuU.U 


a A . 
o 


o.U 




Ft m 1 
[lUlJ 


zuU.U 


Q A . 
o At 


o.U 


PQT 


Fl ni 1 
[lUlJ 


zDU.U 


Si n 

o i^t 


o.U 


PQT 


Fl ni 1 
[lUlJ 


zDU.U 




o.U 


PQT 


Ft ni 1 
[lUlJ 


zDU.U 


S P 

O r 


n 
z.U 


PQT 
1 


Fl ni 1 
[iUiJ 


zoU.U 


Q P 


z.U 




Fl ni 1 
[lUiJ 


zui.UU 


QQ P 
OO 1 


A 1 
4.1 


1 rtiUiVlr 


Fnnl 

[yyj 


261.9 


108 CT 


3.0 


Preiburff PSI 


[100] 


280.0 


109 CT 


3.0 


Freiburg, PSI 


[100] 


300.2 


111 CT 


2.6 


Freiburg, PSI 


[100] 


312.0 


24 P 


4.0 


SATURNE 


[102] 


312.0 


11 


4.0 


SATURNE 


[103] 


318.0 


8 Rt 


3.0 


PSI 


[101] 


318.0 


8 At 


3.0 


PSI 


[101] 


318.0 


5 At 


3.0 


PSI 


[101] 


318.0 


8 Dt 


3.0 


PSI 


[101] 


318.0 


5 Dt 


3.0 


PSI 


[101] 


318.0 


8 P 


2.0 


PSI 


[101] 


318.0 


5 P 


2.0 


PSI 


[101] 


820.1 


110 rr 


2.1 


Froibnrg. PSI 


[100] 


1(1.1) 


112 a 




Fri-ilmrs. I'Sl 


[100] 



"This data set is floated in tlic calculations of table 8 because all current 
phase shift analyses and np potentials predict a norm that is about 4 standard 
deviations oflF the experimental normalization error of 4.9%. 
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Table 10. x'^/datum for various sets of neutron-proton differential cross section 
data below 350 MeV applying some recent phase shift analyses (PSA) and np 
potentials. 





VPI/GWU 


Nijmegen 1993 


Argonne Vis 


CD-Bonn 




PSA [88] 


PSA [64] 


potential [89] 


potential [66] 


Bonner et al. [104], 652 data 


1.18 


1.08 


1.20 


1.10 


Freiburg/PSI [100], 871 data 


7.66 


8.62 


8.58 


8.14 


Uppsala [105, 98], 109 data 


3.40 


6.15 


5.20 


6.41 


Louvain-la-Neuve [96], 84 data 


3.22 


3.15 


3.12 


3.17 



Freiburg data instead, the latter can be reproduced witli a /datum = 2.64 (the 
Bonner data produce /datum = 1.84 for this fit). Thus, the Freiburg data cannot 
be reproduced with the same accuracy as the Bonner data, even if one restricts the 
np a data exclusively to Freiburg. This may be seen as an indication that there are 
inconsistencies within the Freiburg data. Certainly, the Bonner data and the Freiburg 
data are inconsistent with each other. Similar problems arc observed with the Uppsala 
data, which are included in the VPI/GWU analysis and excluded from the Nijmegen 
PSA. The problems with the np differential cross sections deserve further systematic 
investigation. 

6. The new high-precision NN potentials 

In the 1990's, a focus has been on the quantitative aspect of NN potentials. Even the 
best NN models of the 1980's [24, 107] fit the NN data typically with a x^/datum « 2 or 
more. This is still substantially above the perfect x^ /datum w 1. To put microscopic 
nuclear structure theory to a reliable test, one needs a perfect NN potential such that 
discrepancies in the predictions cannot be blamed on a bad fit of the NN data. 

Based upon the Nijmegen analysis and the (pruned) Nijmegen database, new 
charge-dependent NN potentials were constructed in the early/mid 1990's. The groups 
involved and the names of their new creations are, in chronological order: 

• Nijmegen group [108]: Nijm-I, Nijm-II, and Reid93 potentials. 

• Argonne group [89]: Vig potential. 

• Bonn group [65, 66]: CD-Bonn potential. 

All these potentials have in common that they use about 45 parameters and fit the 
(pruned) 1992 Nijmegen data base with a x^/datum « 1. However, as discussed in the 
previous section, since 1992 the pp database has substantially expanded and for the 
current database the x^ /datum produced by some potentials is not so perfect anymore 
(cf. table 8). 

6.1. Theoretical aspects 

Concerning the theoretical basis of these potential, one could say that they are all 
more or less — constructed 'in the spirit of meson theory' (e.g., all potentials include 
the one-pion-exchange contribution). However, there are considerable differences in 
the details leading to considerable off-shell differences among the potentials. 

To explain these details and differences in a systematic way, let us first sketch the 
general scheme for the derivation of a meson-theoretic potential. 
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One starts from field-theoretic Lagrangians for meson-nucleon coupling, which 
are essentially fixed by symmetries. Typical examples for such Langrangians are: 

Cps = - gpsi>il''i>^^'"^ (44) 

Cs = - ff.V^V"/'^'^ (45) 

where ps, s, and v denote pseudoscalar, scalar, and vector couplings/fields, 
respectively. 

The lowest order contributions to the nuclear force from the above Lagrangians 
arc the second-order Fcynman diagrams which, in the center-of-mass frame of the two 
interacting nucleons, produce the amplitude: 

^ ^ ii.(.Ti°'..(c.)P.y(-q;)rr'^,(-q) _ ,„) 

(g ' - qY -m^ 

where F^"^ [i = 1,2) are vertices derived from the above Lagrangians, Ui are Dirac 

spinors representing the nucleons, and q and q' are the nucleon relative momenta in 
the initial and final states, respectively; Pa divided by the denominator is the meson 
propagator. 

The simplest meson-exchange model for the nuclear force is the one-boson- 
exchange (OBE) potential [5] which sums over several second-order diagrams, each 
representing the single exchange of a different boson, a: 

^(q',q) = • (^s) 

a 

As customary, we include form factors, Fa{(^^^)^ applied to the meson-nucleon 
vertices, and a square-root factor M/\/E'E (with E = sj -|- and E' = 
is the nucleon mass). The form factors regularize the amplitudes 
for large momenta (short distances) and account for the extended structure of 
nucleons in a phenomenological way. The square root factors make it possible to 
cast the unitarizing, relativistic, three-dimensional Blankenbecler-Sugar equation for 
the scattering amplitude (a reduced version of the four-dimensional Bcthc-Salpctcr 
equation) into a form which is identical to the (nonrelativistic) Lippmann-Schwinger 
equation (see reference [5] for details). Thus, equation (48) defines a relativistic 
potential which can be consistently applied in conventional, nonrelativistic nuclear 
structure. 

Clearly, the Feynman amplitudes, equation (47), are in general nonlocal 
expressions; i. e., Fourier transforming them into configuration space will yield 
functions of r and r', the relative distances between the two in- and out-going nucleons, 
respectively. The square root factors create additional nonlocality. 

While nonlocality appears quite plausible for heavy vector-meson exchange 
(corresponding to short distances), we have to stress here that even the one-pion- 
exchange (OPE) Feynman amplitude is nonlocal. This is important because the pion 
creates the dominant part of the nuclear tensor force which plays a crucial role in 
nuclear structure. 

Applying F''^^ = 5(^75 in equation (47), yields the Feynman amplitude for one- 
pion exchange, 

. . , , , gl {E' + M){E + M) ( a,.q' a.-q 

*^7r(q ,q) = - 



4M2 (q'-q)2+TO2 \E' + M E + M 
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Figure 10. Half off-shell ^Si-^Di amplitude for the relativistic CD-Bonn 
potential (solid line), equation (48). The dashed curve is obtained when the 

local approximation, equation (51), is used for OPE, and the dotted curve results 
when this approximation is also used for one-p exchange, q' = 265 MeV/c. 



0-2 q' 



, , (49) 

where m.,r denotes the pion mass and isospin factors are suppressed. This is the original 
and correct result for OPE. 

If one now introduces the drastic approximation, 

E' PiEr^iM , (50) 

then one obtains the momentum space representation of the local OPE, 

ll (<Ti.k)(a2-k) 



1^1^°^^ (k) 



4M2 k2 + ml 



V^l^"^^ equals iA^n 



(51) 
Thus, the 



with k = q' — q. Notice that on-shell, i. e., for |q'| — 
nonlocality afTccts the OPE potential off-shell. 

Fourier transform of equation (51) yields the well-known local OPE potential in 
r-space. 



9I / ™^ y 

127r \2m) 



+ 1 



47r 



rriTrr 



(5(3) (r 



-5^12 



(52) 



Notice, however, that this 'well-established' local OPE potential is only an 
approximative representation of the correct OPE Fcynman amplitude. A QED analog 
is the local Coulomb potential versus the full field-theoretic one-photon-exchange 
Feynman amplitude. 

It is now of interest to know by how much the local approximation changes the 
original amplitude. This is demonstrated in figure 10, where the half off-shell ^Si- 
^Di potential, which can be produced only by tensor forces, is shown. The on-shell 
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Table 11. Modern high-precision NN potentials and their predictions for the 
Iwu- and lliree-nucieoii l>ound si ales. 





CD- Bonn 

[66] 


Nijm-I 

[108] 


Nijm-II 
[108] 


Reid93 
[108] 


Vis 
[89] 


Nature 


Character 


nonlocal 


mixed" 


local 


local 


local 


nonlocal 


Deuteron properties: 
Quadr. moment (fm^) 
Asymptotic D/S state 
D-state probab. (%) 


0.270 
0.0256 
4.85 


0.272 
0.0253 
5.66 


0.271 
0.0252 
5.64 


0.270 

0.0251 
5.70 


0.270 

0.0250 
5.76 


0.276(2)'' 
0.0256(4)= 


Triton binding (MeV): 
nonrcl. calculation 
rclativ. calculation 


8.00 
8.2 


7.72 


7.62 


7.63 


7.62 


8.48 



° Central force nonlocal, tensor force local. 

Corrected for meson-exchange currents and relativity. 
= Reference [62]. 



momentum q' is held fixed at 265 MeV/c (equivalent to 150 MeV laboratory energy), 
while the off-shell momentum q runs from zero to 2000 MeV/c. The on-shell point 
{q = 265 McV/c) is marked by a solid dot. The solid curve is the CD-Bonn potential 
which contains the full, nonlocal OPE amplitude equation (49). When the static/local 
approximation, equation (51), is made, the dashed curve is obtained. When this 
approximation is also used for the onc-p exchange, the dotted curve results. It is 
clearly seen that the static/local approximation substantially increases the tensor 
force off-shell. Certainly, we are not dealing here with negligible effects, and the 
local approximation is obviously not a good one. 

Even though the spirit of the new generation of potentials is more sophisticated, 
only the CD-Bonn potential uses the full, original, nonlocal Feynman amplitude for 
OPE, equation (49), while all other potentials still apply the local approximation, 
equations (51) and (52). As a consequence of this, the CD-Bonn potential has a weaker 
tensor force as compared to all other potentials. This is reflected in the predicted 
D-statc probability of the deuteron, Pd, which is due to the nuclear tensor force. 
While CD-Bonn predicts Pd = 4.85%, the other potentials yield Pd = 5.7(1)% (cf. 
table 11). These differences in the strength of the tensor force lead to considerable 
differences in nuclear structure predictions. An indication of this is given in table 11: 
The CD-Bonn potentials predicts 8.00 MeV for the triton binding energy, while the 
local potentials predict only 7.62 MeV. More discussion of this aspect can be found in 
references [5, 109]. 

The OPE contribution to the nuclear force essentially takes care of the long- 
range interaction and the tensor force. In addition to this, all models must describe 
the intermediate and short range interaction, for which very different approaches are 
taken. The CD-Bonn includes (besides the pion) the vector mesons p(769) and u;(783), 
and two scalar-isoscalar bosons, a, using the full, nonlocal Feynman amplitudes, 
equation (47), for their exchanges. Thus, all components of the CD-Bonn arc nonlocal 
and the off-shell behavior is the original one that is determined from relativistic field 
theory. 

The models Nijm-I and Nijm-II are based upon the Nijmegen78 potential [81] 
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Tiab T,^^ (MeV) 

Figure 11. Upper part: Matrix elements V{q',q) of tlie ^Sq and potentials 
for the CD-Bonn (solid line), Nijm-I (dashed), Nijm-II (dash-dot), Argonne Vis 
(dash-triple-dot) and Reid93 (dotted) potentials. The diagonal matrix elements 
with q' = q = 265 MeV/c (equivalent to Ti^f, = 150 MeV) are marked by a solid 
dot. The corresponding matrix element of the scattering ii"-matrix is marked by 
the star. Lower part: Predictions for the np phase shifts in the ^So and ^Si 
state by the five potentials. The five curves are essentially indistinguishable. The 
solid dots represent the Nijmegen multi-energy np analysis [64]. 

which is constructed from approximate OBE ampUtudes. Whereas the Nijm-II uses 
the totally local approximations for all OBE contributions, the Nijm-I keeps some 
nonlocal terms in the central force component (but the Nijm-I tensor force is totally 
local). Nonlocalities in the central force have only a very moderate impact on nuclear 
structure as compared to nonlocalities in the tensor force. Thus, if for some reason 
one wants to keep only some of the original nonlocalities in the nuclear force and not 
all of them, then it would be more important to keep the tensor force nonlocalities. 

The Rcid93 [108] and Argonne Vis [89] potentials do not use meson-exchange for 
intermediate and short range; instead, a plicnomenological parametrization is chosen. 
The Argonne Vis uses local functions of Woods-Saxon type, while Reid93 applies 
local Yukawas of multiples of the pion mass, similar to the original Reid potential of 
1968 [110]. At very short distances, the potentials are regularized either by exponential 
{Vi8, Nijm-I, Nijm-II) or by dipole (Reid93) form factors (which are all local functions). 
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In figure 11, the five high-precision potentials (in momentum space) and their 
phase shift predictions arc shown, for the ^S^ and states. While the phase shift 
predictions are indistinguishable, the potentials differ widely — due to the theoretical 
and mathematical differences discussed. Note that NN potentials differ the most in 
S'-wavcs and converge with increasing L (where L denotes the total orbital angular 
momentum of the two-nucleon system). 

6.2. Charge dependence 

All new potentials are charge-dependent which is essential for obtaining a good for 
the pp and up data. Thus, each potential comes in three variants: pp, np, and nn. 
All potentials include the CIB effect from OPE. However, as discussed in section 

2.2, pion mass splitting creates further CIB effects through the diagrams of 27r 
exchange and other two-boson exchange diagrams that involve pions. Another source 
of CIB is irreducible exchange. Recently, these contributions have been evaluated 
in the framework of chiral perturbation theory by van Kolck et al. [34] . In L > states, 
the size of this contribution is typically the same as the CIB effect from TBE. Thus, 
TBE and 7r7 create sizable CIB effects in states with X > 0. Therefore, a thoroughly 
constructed, modern, charge-dependent NN potential should include them. The NN 
potentials [108, 89] ignore these contributions while the latest CD-Bonn update [66] 
incorporates them. 

A similar comment can be made about CSB. Most potentials include; only the 
most trivial effects from nucleon mass splitting, namely the effect on the kinetic energy 
and on the OBE diagrams. However, as discussed in section 2.1.2, there are relatively 
large contributions from TBE that fully explain the CSB scattering length difference. 
Because of the outstanding importance of the CSB effect from TBE, it should be 
included in NN force models (and, therefore, it has been incorporated in the latest 
update of the CD-Bonn potential [66]). To have distinct pp and nn potentials is 
important for addressing several interesting issues in nuclear physics, like, the '^H-^He 
binding energy difference and the Nolen-Schiffer (NS) anomaly [111] regarding the 
energies of neighboring mirror nuclei. Potentials that do not include any CSB have 
no chance to ever explain these phenomena. Some potentials that include CSB focus 
on the ^Sq state only, since this is where the most reliable empirical information is. 
However, even this is not good enough. A recent study [112] has shown that the 
CSB in partial waves with L > as derived from the Bonn model is crucial for a 
quantitative explanation of the NS anomaly. 

6.3. Extrapolating low-energy potentials towards higher energies 

NN potentials designed for nuclear structure purposes are typically fitted to the NN 
scattering data up to pion production threshold or slightly beyond (e. g., 350 MeV). 
A very basic reason for this is that a real potential cannot describe the inelasticities of 
particle production. On the other hand, nuclear structure calculations are probably 
sensitive to the properties of a potential above 350 MeV. For example, the Brueckner 
G-matrix, which is a crucial quantity in many microscopic approaches to nuclear 
structure, is the solution of the integral equation. 




(53) 
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Figure 12. The £2 mixing parameter and the ^F:i phase shift up to 1000 MeV 
lab. energy for various potentials as denoted (N-I and N-II refer to the Nijmegen 
potentials). Solid dots represent the Nijmegen PSA [64] and open circles the 
VPI/GWU analysis SM99 [72]. 



(where M* denotes the effective nucleon mass and Q the Pauh projector). Notice 
that the potential V is involved in this equation for all momenta from zero to infinity, 
on- and off-shell. Now, it may very well be true that, as the momenta increase, their 
importance decreases (due to the short-range repulsion of the nuclear force and the 
associated short-range suppression of the nuclear wave function). However, it is also 
true that the impact of the potential does not suddenly drop to zero as soon as the 
momenta involved become larger than the equivalent of 350 MeV lab. energy. Thus, 
there are good arguments why NN potentials should extrapolate in a reasonable way 
towards higher energies. 

We have investigated this issue and found good and bad news. The good news 
is that most potentials reproduce in most partial waves the NN phase shifts up to 
about 1000 MeV amazingly well. The bad news is that there are some singular cases 
where the reproduction of phase parameters for higher energies is disturbingly bad. 
The two most notorious cases are shown in figure 12. Above 350 MeV, the £2 mixing 
parameter is substantially underpredicted by both Nijmegen potentials (N-I and N-II). 
The reason for this is that, for 62, both potentials follow very closely the Nijmegen 
PSA [64] (solid dots in figure 12) up to 350 MeV. Thus, these potentials are faithfull 
extrapolations of the Nijmegen PSA to higher energies. Since this extrapolation is 
wrong, the suspicion is that the Nijmegen PSA has a wrong trend in the energy range 
250-350 MeV. New data on pp spin transfer coefficients [61] in the energy range 300-500 
MeV could resolve the issue. 

A similar problem occm's in ^F^ (figure 12). Here, the dashed curve (N-I) is the 
extrapolation of the Nijmegen PSA, indicating that the analysis may have the wrong 
trend in the energy range 200-350 MeV. 
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We note that, in the two channels discussed, the inelasticity has little impact on 
the phase parameters shown and would not fix the problems. 

The moral is that one should not follow just one analysis, particularly, if that 
analysis is severely limited in its energy range. It is important to also keep the broad 
picture in mind. 

7. The theory of nuclear forces: future directions 

7.1. Critical summary of current status 

During the past decade or so, the research on the NN interaction has proceeded 
essentially along two lines. There was the phenomenological line which has produced 
the high-precision, charge-dependent NN potentials [108, 89, 65, 66]. This was 
practical work, necessary to provide reliable input for exact few-body calculations 
and nuclear many-body theory. 

The goal of the second line of research was to approach the problem on a more 
fundamental level. Since about 1980, we have seen many efforts to derive the nuclear 
force from the underlying theory of strong interactions, quantum chromo dynamics 
(QCD). Due to its nonperturbative character in the low-energy regime, QCD cannot be 
solved exactly for the problem under consideration. Therefore, so-called QCD-related 
or QCD- inspired models have been developed, like, Skyrmion or Soliton models [113] 
and constituent quark cluster models. Among the quark models, one may distinguish 
between two types: the hybrid models [114] that include meson exchange and the 
quark delocalization and color screening models [115] that do not need (and do not 
include) meson exchange to create the intermediate-range attraction of the nuclear 
force. The success of all these efforts based upon QCD modelling is mixed. The 
qualitative features of the nuclear force are, in general, predicted correctly, but none 
of the models is sufficiently quantitative such that it would make sense to apply it in 
miscroscopic nuclear structure calculations. 

In summary, one problem of the current status in the field is that quantitative 
models for the nuclear force have only a poor theoretical background, while theory 
based models yield only poor results. This discrepancy between theory and practice 
has become rather larger than smaller, in the course of the 1990s. Another problem is 
that the 'theory based models' are not strictly derived from QCD, they are modeled 
after QCD — often with handwoven arguments. Thus, one may argue that these 
models are not any better than the traditional meson-exchange models (that are 
nowadays perceived as phenomenology). The purpose of physics is to explain nature 
in fundamental terms. The two trends just discussed are moving us away from this 
aim, which is reason for serious concern. 

Therefore, the main goal of future research on the nuclear force must be to 
overcome the above discrepancies. To achieve this goal, we need a basic theory that is 
amenable to calculation and yields quantitative results. 

7.2. The effective field theory concept 

In recent years, the concept of effective field theories (EFT) has drawn considerable 
attention in particle and nuclear physics [116, 117, 118]. The notion of effective field 
theories may suggest a difference to fundam,6ntal field theories. However, it is quite 
likely that all field theories (including those that we perceive presently as fundamental) 
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are effective in the sense that they are low-energy approximations to some 'higher' 
theory. 

The basis of the EFT concept is the recognition of different energy scales in 
nature. Each energy level has its characteristic degrees of freedom. As the energy 
increases and smaller distances arc probed, new degrees of freedom become relevant 
and must be included. Conversely, when the energy drops, some degrees of freedom 
become irrelevant and are frozen out. 

To model the low-energy theory, one relies on a famous 'folk theorem' by 
Weinberg [119, 120] which states: 

If one writes down the most general possible Langrangian, including all terms 
consistent with assumed symmetry principles, and then calculates matrix 
elements with this Langrangian to any given order of perturbation theory, 
the result will simply be the most general possible S-matrix consistent with 
analyticity, perturbative unitarity, cluster decomposition, and the assumed 
symmetry principles. 

The essential point of an effective field theory is that we are not allowed to make any 
assumption of simplicity about the Lagrangian and, consequently, we are not allowed 
to assume renormalizability. The Langrangian must include all possible terms, because 
this completeness guarantees that the effective theory is indeed the low-energy limit 
of the underlying theory. Now, this implies that we are faced with an infinite set 
of interactions. To make the theory managable, we need to organize a perturbation 
expansion. Then, up to a certain order in this expansion, the number of terms that 
contribute is finite and the theory will yield a well-defined result. 

In strong interactions, the transition from the 'fundamental' to the effective level 
happens through a phase transition that takes place around Aqcd ~ 1 GeV via 
the spontaneous breaking of chiral symmetry which generates pscudoscalar Goldstone 
bosons. Therefore, at low energies {E < Aqcd), the relevant degrees of freedom 
are not quarks and gluons, but pseudoscalar mesons and other hadrons. Approximate 
chiral symmetry is reflected in the smallness of the masses of the pseudoscalar mesons. 
The effective theory that describes this scenario is known as chiral perturbation theory 
(XPT) [117, 121, 122]. 

If we believe in the basic ideas of EFT, then, at low energies, xPT is as 
fundamental as QCD at high energies. Moreover, due to its perturbative arrangement, 
xPT can be calculated: order by order. So. here we may have what we are asking for: 
a basic theory that is amenable to calculation. Therefore, xPT has the potential to 
overcome the discrepancy between theory and practice that has beset the theoretical 
research on the nuclear force for so many years. 

7.3. Chiral perturbation theory and nuclear forces 

The idea to derive nuclear forces from chiral effective Lagrangians is not new. 
A program of this kind was started some 10 years ago by Weinberg [123, 124], 
Ordoiiez [125], and van Kolck [126, 127, 128] which produced the first chiral NN 
potential by Ordonez, Ray, and van Kolck [129]. 

After the program was initiated, considerable activity ensued [130, 131, 132, 133, 
134, 135, 136, 137, 138, 139, 140, 141]. For a recent review, see Ref. [142]. Even though 
all authors start from chiral effective Langrangians, there are differences in the details. 
There is, for example, the KSW scheme [135] in which the amplitude of interest is 
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calculated perturbatively. On the other hand, Weinberg proposed to use xPT ^ot 

computing the NN potential which consists of irreducible diagrams. The S'-matrix 
is then generated by the Schrodinger equation. The first comprehensive work using 
the Weinberg scheme was done by Ordonez, Ray, and van Kolck [f 25, 126, 129] who 
applied time-ordered ('old-fashioned') perturbation theory to calculate the irreducible 
diagrams that define the potential. This potential posesses one unpleasant property, 
namely, it is energy dependent. One can avoid this problem by using the method 
of unitary transformations, a method that was pioneered by Okubo [143]. The 
Okubo tranformation is applied in the recent work by Epelbaum, Glockle, and 
Meifiner [144, 145] who construct chiral NN potentials in leading order (LO) of 
XPT, next-to-leading order (NLO), and next-to-next-to-leading order (NNLO). A 
systematic improvement in the ability of the model to reproduce the NN data is 
obseerved when stepping up the orders of the chiral expansion. The NNLO potential 
of Epelbaum et al. [145] describes the np phase shifts well up to about 100 McV; 
above this energy there are discrepancies in some partial waves. This most recent 
chiral NN potential represents great progress as compared to earlier ones, however, 
for meaningfuU applications in microscopic nuclear structure, futher quantitative 
improvements are necessary. 

There is one particularly attractive aspect to the xPT approach in regard to 
those nuclear structure applications. If, in the traditional approach, one wants to 
reproduce, e. g., the experimental binding energies of the triton, the alpha particle 
or other nuclei, one complements the NN potential with a (phenomenological) three- 
nuclcon force (3NF) [146]. Since different NN potentials leave different discrepancies to 
experiment (cf. table 11), the 3NF is adjusted from potential to potential. From a more 
fundamental point of view, this proceedure is very unsatisfactory, since it lacks any 
underlying systcmatics. However, within the framework of traditional meson theory, 
there is nothing else you can do, because there is no a priori connection between the 
off-shell NN potential and the existence of certain many-body forces. 

In the framework of xPT. there is this connection from the outset. In each order 
of xPT, the two-nucleon force is well-defined on- and off-shell and it is also well-defined 
which 3NF terms occur in that order. At least that's how it should work 'in theory'. 
How it works out in practise remains to be seen. 

8. NN scattering at intermediate and high energies 

Even though the focus of this review is on the NN interaction at low energy (below 
pion production threshold), we like to give an indication of the exciting developments 
at higher energies to ensure that the reader gets an impression of the broader picture. 

In the 1990's, a new experimental effort to measure pp scattering observables 
in the range 0.5 to 2.5 GeV of projectile energy was started. The experiments are 
conducted by the EDDA group [147, 148] at the cooled proton synchrotron (COSY) 
at Jiilich, Germany, and use an internal target. The outstanding features of this new 
generation of experiments are high precision (high statistics) and small energy steps. 
Other experiments on NN observables at intermediate and high energies have been 
performed at SATURN II (Saclay, France) [149] and at facilities around the world; 
for a complete listing, see Ref. [150]. Moreover, there are the pp analyzing powers 
and spin correlation parameters, Cnn, that were measured up to lab. energies of 
11.75 GeV by Alan Krisch and coworkers some 20 years ago at the Argonne Zero- 
Gradient-Synchrotron (ZGS) [151]. The latter data were never explained in a really 
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satisfatory way in spite of some efforts [152, 153]. During the past decade, theory 
has shown Uttlc interest in NN scattering at energies of 1-10 GcV and only a few 
papers can be cited [154, 155, 156, 157, 158, 159]. The reason for the neglect is 
probably that these energies are too high for traditional nuclear physicists and too 
low for high energy physicists; so, this range is the stepchild of both professions. But 
this is what makes this window particularly interesting: the energies are too high 
for xPT and too low for perturbative QCD. This fact implies that we have to find an 
appropriate phenomenology which calls for some phantasy and creativity. Relativistic, 
chiral meson models that include heavy mesons [like, p(770) and a;(782)] and nucleon 
resonances [like, the A(1232) isobar] are an obvious choice. In the language of EFT, 
these models can be justified on the basis of 'resonance saturation'. It will also be 
interesting to attempt matching of these high energy models with the xPT model 
discussed in the previous section. 

It is well-known that relativistic meson models work satisfactorily up to about 
1.2 GeV [160, 161, 162, 163, 164]; for a summary and critical discussion see section 7 
of Ref. [5]. However, at energies above 1.5 GeV, characteristic problems occur, some 
of which are [165]: 

• The predicted elastic NN cross sections are too large and grow with energy while 

experimentally they drop. 

• The pp analyzing powers are predicted too large and for other spin observables 
(like, Cjviv) even the sign is predicted wrong. 

The first problem listed above is well-known since the late 1950's. The amplitude 
produced by vector-meson exchange is proportional to s/ (f — m^) (with rUy the vector- 
meson mass and s, t the usual Mandelstam variables) which causes the elastic cross 
sections to rise with energy. Because of this problem, Regge theory [166, 167, 168, 169] 
was invented. Concerning the second problem listed above, we do not have a clue at 
this time. When the vector-meson contribution is 'phased out' such as to produce the 
correct elastic cross sections, then the problem with the analyzing powers persists, 
which is not what we expected. In any case, this energy region offers a wealth of good 
and unexplained data and a great diversity of potentially appropriate models, since 
we are truely at the intersection of nuclear and particle physics. 

9. Summary and outlook 

In the 1990's, we made major progress in our grasp on the nuclear force. NN data 
of unparalleled precision were produced, particularly at TUNL [90, 91], lUCF [61, 
84, 85, 86, 87], and COSY [147, 148]. The art of NN phase shift analysis advanced 
substantially [64]. Based upon this empirical progress, NN potentials of unprecedented 
accuracy (x^/datum « 1) were contructed [108, 89, 65, 66]. As a consequence of this, 
exact few-body calculations and miscroscopic; nuclear many-body theory can now be 
based upon input that is more reliable than ever. 

In spite of all this good news, there are also several questions concerning the NN 
interaction that remain open and future research should continue to address them. 
Among the more technical problems is the neutron-neutron scattering length where 
experiments are still in contradiction. There are problems with the determination 
of a precise value for the ttNN coupling constant and, not unrelated, the huge and 
expensive database of neutron-proton differential cross sections [104, 105, 84, 100] is 
inconsistent. 
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On the theoretical side, we are still lacking a derivation of the nuclear force that 
is based upon theory (in the true sense of the word) and produces a quantitative NN 
potential. Moreover, our understanding of charge-dependence of the NN interaction 
is still incomplete, since we are not able to explain 25% of the charge-dependence of 
the ^So scattering length. 
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